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Stationary and axisymmetric ideal magnetohydrodynamic (MHD) accretion onto a black hole is 
studied analytically. The accreting plasma ejected from a plasma source with low velocity must 
be super-fast magnetosonic before passing through the event horizon. We work out and apply a 
trans-fast magnetosonic solution without the detailed analysis of the regularity conditions at the 
magnetosonic point, by introducing the bending angle j3 of magnetic field line, which is the ratio of 
the toroidal and poloidal components of the magnetic field. To accrete onto a black hole, the trans- 
magnetosonic solution has some restrictions on (3, which are related to the field-aligned parameters 
of the MHD flows. One of the restrictions gives the boundary condition at the event horizon for 
the inclination of a magnetic field line. We find that this inclination is related to the energy and 
angular momentum transport to the black hole. Then, we discuss the spin- up/down process of a 
rotating black hole by cold MHD inflows in a secular evolution timescale. There are two asymptotic 
states for the spin evolution. One is that the angular velocity of the black hole approaches to that of 
the magnetic field line, and the other is that the spin-up effect by the positive angular momentum 
influx and the spin-down effect by the energy influx (as the mass-energy influx) are canceled. We 
also show that the MHD inflows prevents the evolution to the maximally rotating black hole. 



PACS numbers: 97.60.Lf, 95.30.Qd, 95.30.Sf 



I. INTRODUCTION 



As central engines of active galactic nuclei, some com- 
pact X-ray sources and gamma-ray bursts, black hole ac- 
cretion systems are generally accepted. Recently, it is 
also confidently expected that the magnetic field is con- 
sidered to play a very important role in the black hole ac- 
cretion systems |7| . In addition, Meier [2j proposed that 
the inner part of a black hole accretion inflow may enter 
a magnetically-dominated phase. Although black holes 
cannot support a magnetic field by themselves, they can 
be immersed in the magnetic field generated by currents 
in black hole accretion flows. Then, above an equatorial 
disk and between the disk and a black hole, we can expect 
a magnetically-dominated "magnetosphere" like our Sun 
or pulsars, where large-scale (global) well-ordered mag- 
netic fields extend to distant regions, and some part of 
the fields extend inwardly and enter the black hole. In 
this magnetosphere, more magnetically important phe- 
nomena are expected. 

The research of most pulsating X-ray sources as accret- 
ing neutron stars is based on the qualitative features of 
their spectra and their spin- up/down rates. To discuss 
these period changes in these sources, Ghosh & Lamb 
[H studied the accretion torque on a magnetic neutron 
star from a Keplerian disk by solving the hydromagnetic 
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equations along the dipole magnetic field lines. In the 
black hole case, unfortunately, the spin of a black hole is 
unclear for every black hole candidate. This is because 
many uncertain factors about the black hole-accretion 
disk systems still remain to be understood; e.g, the black 
hole mass (the size of the black hole) , the distribution of 
magnetic field and accreting plasma, and the radiative 
process by the accreting plasma, etc. The spin of a black 
hole would be of basic importance for considering the 
environment around the black hole and the activities of 
the black hole-accretion disk systems. Then, the under- 
standing of the process of the hole's spin change would 
provide an important tool for exploring the astrophysical 
activity of a black hole accretion system. 

The evolution of black hole spin by accreting matter 
was first studied by Bardeen Q. Assuming that the gas 
trickles down the event horizon from the inner-edge of an 
equatorial accretion disk, which was set to be located at 
the last stable circular orbit of the Kerr geometry, it was 
shown that a non-rotating black hole would always be 
span up during accretion. However, Thorne [eJ showed 
that photon capture by a gas accreting black hole would 
become important in the last stages of spin-up, where 
photons are emitted by the accretion disk, and reported 
a limiting value of a/m ~ 0.9982, where m and a de- 
note the mass and angular momentum per unit mass of 
the black hole, respectively. This is because the hole's 
capture cross-section is greater for photons of negative 
angular momentum than photons of positive angular mo- 
mentum. Then, the captured photons prevent complete 
spin-up to a/m = 1. When we consider a black hole 



2 



immersed in a magnetosphere, the energy and angular 
momentum extraction due to a magnetic torque on the 
hole can be expected [f| 0] ■ The magnetic field also pre- 
vents the spin-up to a/m = 1. 

The MHD inflow carries the total energy E and angular 
momentum L into a black hole. Then, the black hole will 
obtain mass 8m = E and angular momentum 5 J = L, 
where J = am. Although in this paper we assume that 
the background metric does not change by the accretion 
process (at least within our observational timescale), the 
black hole will evolve in a secular timescale. That is, for 
example, the angular velocity, the spin and the surface 
area will be modified. From the hole's angular velocity 
loh = a/2mrH, where ru = m + \J m 2 — a 2 is the hori- 
zon's radius, we can obtain the relation 
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H 



2mr jjV 'm 2 — a 2 



L 



-(m + r H )E 



(1) 



We may expect that the magnetic field lines connected to 
the lower latitude region of the event horizon contribute 
to the hole's spin-up and the field lines connected to the 
higher latitude (polar) region of the event horizon con- 
tribute the spin-down (e.g., see H). To discuss the spin 
evolution, we need to construct a reasonable model of 
MHD accretion in a black hole magnetosphere. 

In this paper, we consider a stationary and axisymmet- 
ric black hole magnetosphere with a magnetized accre- 
tion disk, where plasma flows stream along magnetic field 
lines under the ideal MHD approximation. The global 
magnetic fields can be generated by current in the equa- 
torial magnetized accretion disk and its off-equatorial 
plasma in the disk's corona region. In this black hole 
magnetosphere, the magnetic field lines originated from 
the inner part of the accretion disk (within several grav- 
itational radius) connect to the event horizon, and the 
fields from the outer part of the accretion disk extend to 
distant regions. The plasma ejected from the disk surface 
must stream inwardly toward the event horizon along the 
black hole-disk magnetic field lines or outwardly toward 
distant regions along the disk's open magnetic field lines; 
so that, we need both MHD "inflow" and "outflow" so- 
lutions, which are obtained from the general relativistic 
Bernoulli equation. 

The general relativistic Bernoulli equation for magne- 
tized flows, which is the relativistic extension from the 
work by Weber & Davis @, was formulated by Camen- 
zind and Takahashi et al.Q- We consider physically 
acceptable MHD accretion solutions that start from a 
plasma source with low velocity, which is sub-Alfvenic ve- 
locity (or sub-slow magnetosonic when we consider a hot 
flow). The physically acceptable accelerated inflow falls 
into the black hole after passing through the Alfvcn point 
and the fast magnetosonic point. At the magnetosonic 
points, which are X-type (saddle type) critical points, 
a smooth transition from sub-magnetosonic to super- 
magnetosonic flows must occur. Takahashi inves- 
tigated the critical conditions on the trans-magnetosonic 
flows in the Kerr geometry, and presented the restriction 



on the field-aligned flow parameters for the appearance 
of the magnetosonic point on the accretion solutions. In 
general, the task for finding a trans-magnetosonic solu- 
tion is troublesome because the relations between the 
flow parameters are very complicated. Therefore, in this 
paper, we will study the trans-fast magnetosonic MHD 
accretion solution without the detailed analysis of the 
regularity conditions at the magnetosonic point, by in- 
troducing the bending angle j3 of magnetic field, which 
is defined as the ratio of the toroidal and poloidal com- 
ponents of the magnetic field. The singular term in the 
equation of motion of MHD plasma is incorporated in 
this function. We only assume that the function (3 is reg- 
ular at the magnetosonic point, and then we find that 
the regularity condition of the trans-magnetosonic flow 
solution is automatically satisfied. Thus, we can easily 
find a trans-magnetosonic solution, when the function 
(3(r,9) is given as a model. In former treatments about 
trans-magnetosonic solution (including the critical point 
analysis), the toroidal component of magnetic field is ob- 
tained by solving the relativistic Bernoulli equation under 
a given poloidal magnetic field. Although, in this paper, 
we assume a function (3 and discuss the restrictions on 
it by field-aligned flow parameters (see Appendix [XJ , we 
can say that our approach is basically equivalent to the 
discussion about the toroidal magnetic field along the 
trans-magnetosonic flow in the former treatments. 

Although we propose a new approach to solve the 
trans-magnetosonic solutions without the complicated 
regularity conditions at the magnetosonic points, the task 
for finding the function — f3(r, 9) is another difficult 
problem, which is not discussed here. In this paper, a 
flux tube is considered, but the cross-section of the mag- 
netic flux tube is not specified beforehand. Of course, 
we should give a consistent function (3(r, 9) along a mag- 
netic field line for a black hole magnetosphere considered. 
To do this, we need to treat the force-balance between 
magnetic field lines streaming the trans-magnetosonic 
flows. However, by considering the relativistic Bernoulli 
equation, we can discuss the restrictions on the trans- 
magnetosonic flows in a black hole magnetosphere with- 
out a complete field configuration satisfying the force- 
balance equation (the relativistic Grad-Shafranov (GS) 
equation [12j|). The GS equation would give another re- 
strictions on flow solutions in addition to the restrictions 
discussed in this paper. This method for black hole accre- 
tion is the general relativistic extension from the special 
relativistic outgoing trans-magnetosonic flows by Tomi- 
matsu & Takahashi [l3|. 

In § HI we review the basic equations for trans- 
magnetosonic flows in Kerr geometry (see also @, Hoj ) 
and introduce the function [3. By assuming a smooth 
function of j3, we show trans-magnetosonic flows without 
the regularity condition analysis. In § IIIII we solve cold 
trans-fast magnetosonic flow solutions. To fall into the 
black hole, the trans-fast magnetosonic solution must sat- 
isfy some conditions for the field-aligned flow parameters. 
We show these restrictions (the necessary condition) in 
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§ IIVI Next, in § [V] we present the boundary conditions 
at the event horizon for ingoing MHD flows. We show 
that the magnetic field configuration restricts the signa- 
ture of energy and angular momentum fluxes across the 
event horizon. To specify the energy and angular mo- 
mentum of MHD inflows, a model of the plasma source 
is necessary. But we can discuss the spin- up/down (as 
a secular evolution) of a rotating black hole by consider- 
ing the restriction on MHD inflows at the event horizon. 
Finally, we give brief remarks in § IVII 



II. BASIC EQUATIONS FOR MHD FLOWS 

We consider stationary and axisymmetric ideal MHD 
flows in Kerr geometry. The background metric is written 
by the Boyer-Lindquist coordinates with c = G = 1, 



ds 2 = 



2mr\ Aamr sin 2 9 
1 — I dtf H dtd(f> 



A sin 2 



^dr 2 
A 



Y,d9 2 



(2) 



where A = 

J2 i „2\2 



r 2 — 2mr + a 2 , X = r 2 + a 2 cos 2 9, A = 
(r 2 + a 2 ) 2 — Aa 2 sin 2 9. The particle number conser- 
vation is (nu a ) :a = 0, where n is the number den- 
sity of the plasma and u a is the fluid 4-velocity. The 
ideal MHD condition is u^F a p = 0, where F pv is the 
electromagnetic tensor. The relativistic Polytropic rela- 
tion is P = KpQ, where po — nm pavt is the rest mass 
density, m pa rt is the mass of the plasma particle and 
r is the adiabatic index. The equation of motion is 
T a ^ : p = 0. The energy- momentum tensor is given by 



rpa/3 



TlflU U 



■ Pg al3 + (1/4tt) [F a x F x/3 + (l/Tjg^F 2 ] , 



where p = (p + P)/n is the relativistic enthalpy and 
F 2 = F^ v F pv . The magnetic and electric fields seen by 
a distant observer are defined by B a = (l/2)rj a i3 1 sk l3 F~ fS 
and E a = F a pk", where k a = (1,0,0,0) is the timelike 
Killing vector and r] a /s 7 s = (—gY^^ap^s- The poloidal 
component B p of the magnetic field seen by a lab-frame 
observer is defined by 



B 2 = -B A B A /G 2 (A = r,0) 

= -\r{d^) 2 + 9 s \d e m/p 
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w ? 



(3) 
(4) 



where ^f(r, 9) is the magnetic stream function (the <fr com- 
ponent of the vector potential, A$). The poloidal com- 
ponent Up of the velocity is defined by u 2 = —u a ua- \ 
Here, we set u p > for both ingoing flows (u r < 0) and 
outgoing flows (u r > 0). ] 

The ideal MHD flows stream along the magnetic field 
lines (i.e., ^(r, 9) — constant lines), and have five field- 
aligned flow parameters; that is, the angular velocity 
of field lines rV('i') = —Ft a/ F^a, the number flux 
per magnetic flux f]^) — nu p /B p , the total energy 
E(^f) = fiut — QfBa,/ (Airrf), the total angular momen- 
tum L(^B) = —puA, — Ba,/ (Atttj) and the entropy, which is 



related to K(^). The relativistic Alfven Mach number 
M is defined by 



M 2 = ^^U 2 p _ pup 



B 1 
p 



B„ 



(5) 



where the term B p = B p /(4irp c ri) is introduced to nondi- 
mensionalize. The relativistic enthalpy can be expressed 



in terms of u p and B p 
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= JL = i+ (E™L) (El 
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1 + Phot 



r-i 



(6) 

where p c = m part and /i ho t = {p- m ]l ^^{^^cvf 1 . The 
term p- ln j is evaluated at the plasma injection point by 
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(7) 

The toroidal magnetic field Ba, = (Asm9/T,)Fe r can be 
expressed in terms of flow's parameters and the Alfven 
Mach number, and the non-dimensional toroidal mag- 
netic field Ba, is defined as 



B, 



1 



B,i 



G<j,E + GtL 
p w ) A-npcT] p w (M 2 - a) 



(8) 



where E = E/p c , L = L/p c , a = g tt + 2g t A,Q F + g^ttp, 
G<t> = 9t<t> + 94>4>^f, and G t = gtt + 9t<t>^F- The locations 
of a(r; ty) = give the inner and outer light surfaces for 
magnetic field lines of flp — O ^ (S&) . When M 2 = a at 
some location, it seems that the toroidal magnetic field 
diverges. Such location is called the Alfven point, be- 
cause the poloidal velocity of the flow equals the Alfven 
wave speed there. To realize the physical trans- Alfvenic 
flow, we also require the condition L/E = —(G^/G^a 
there (see 0]), where the label "A" indicates quantities 
at the Alfven radius. When we set B p > (i.e., de^ > 0) 
in the Northern hemisphere and B p < in the Southern 
hemisphere, we obtain rj > and rj < in the respective 
hemispheres. Then, the direction of the toroidal mag- 
netic field Ba, is also reversed for the equatorial plane. 

The poloidal equation (the relativistic Bernoulli equa- 
tion) that gives the Mach number of the streaming 
plasma in the magnetosphere can be written as 



i 2 



p 2 a- M 4 (aB 2 p 



Bl) 







(9) 



where e = E — flpL. The differential form of the poloidal 
equation ([9]) can be written as 



(In Up)' 



where H 



In 



B, 



(2/p 2 ) (e 2 



A/ 4 



p 2 a 



B, 



CLe 2 ) 



1 



H 



) i is the sound four- velocity and a. 
pT 1 )] 1 / 2 is the sound three-velocity. 



A 

(10) 

iW = [(r_-i)(i- 

The prime is a 



4 



derivative along a stream line d r + (B e / B r )de . It seems 
that the singularities at the Alfven point and the fast and 
slow magnetosonic points are removed from the differen- 
tial form of the poloidal equation, but we should note 
that the critical behavior is included in the term of (Si)'. 
In fact, by substituting expressed as equation (jHJ for 
equation (|10p , we obtain the traditional expression; that 
is, (In tip)' = Af/T>, where the numerator M and the de- 
nominator V cx (up — uaw) 2 (u p — ufm)(u p — usm) are the 
functions of M 2 , r and ^ with field- aligned flow parame- 
ters (see 0] and [ll[ for the analysis at the magnetosonic 
points). The terms waWi «fm, "sm are the Alfven wave 
speed, the fast and slow magnetosonic wave speeds, re- 
spectively. 

When we assume a magnetic flux function *B = if?(r, 9) 
and try to solve the poloidal equation, we need to specify 
a set of five field-aligned parameters satisfying the criti- 
cal conditions at the X-type magnetosonic points and the 
Alfven point. To obtain a physical trans-magnetosonic 
flow solution, the fine-tuning of the parameters is re- 
quired. In general this task is very complicated. How- 
ever, we now propose a new analytical method to study 
the trans-magnetosonic flows without the critical condi- 
tions. We can relate the toroidal magnetic field to the 
poloidal magnetic field by defining the bending angle of 
a magnetic field line as 



P(r,9) 



(11) 



Then, in the differential form of the poloidal equa- 
tion ([101), (P 2 Y takes the place of (S 2 ,)'. When (3 is a 
smooth function at the magnetosonic points, which is 
a natural situation on accretion problems, there is no 
need to analyze the regularity condition there. That 
is, in the r-u p plane, the inclination (u p )' of a flow so- 
lution is determined with a finite value anywhere. At 
the event horizon, the toroidal magnetic field becomes 

Bf = {-g^/^iuJH - n F ){d e *)H and *(r ff , 6) = fi- 
nite, which are the boundary conditions there. The label 
"H" indicates quantities at the event horizon. Then, we 

1 /2 

find the condition (3h = (~.90</>)# y^H — &f)- Thus, we 
can obtain physical accretion solutions passing through 
the event horizon with a finite Mach number. 

We also introduce the poloidal electric-to-toroidal 
magnetic field amplitude ratio £ 2 seen by a zero angu- 
lar momentum observer (ZAMO) as (see also [Hj]) 



E p 



3 

r2 



Bp 
Bt 



(12) 
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where we use the following relations E p = —E Ea = 
[G<f>/ Pw) 2 B p , Bp = —B Ba — ot 2 z Bp, B T = —B^B^ = 
(B^/p w ) 2 . The magnetic and electric fields seen by a 
ZAMO are defined by B a = (l/2)ri a(3lS h' 3 F~< 5 and E a = 
F^h 13 , where = (h%, 0, 0, h%) = a^^.O.O^) is the 



[3). The term a z = \/{g u ) 1/2 = (Y.A/A) 1 / 2 is the 
lapse function and u> = —gt^jg^ is the angular velocity 
of a ZAMO with respect to a distant observer. Then, we 
obtain the relation £ 2 = —g^^p — uj) 2 / j3 2 . At the event 
horizon, we obtain £jy = 1. Although the definitions and 
formalisms introduced in this section are available for hot 
MHD flows (P ^ 0), in the following section for the sake 
of simplicity let's consider the cold MHD flows (P = 0). 



III. COLD TRANS-FAST MAGNETOSONIC 
FLOW SOLUTIONS 

For a cold MHD flow (/} = 1) the poloidal equation © 
is reduced to 



AM - 2BM 2 + C = , 



where 
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■{G^E + GtLf 
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(13) 



(14) 

(15) 
(16) 



with z = -(k + l)pl/(G<t,E + G t L) 2 , and k = (g^E 2 + 
2g t( j ) EL + g u L 2 )/ ' p 2 w . Then, the Alfven Mach number of 
the flow is solved by 



Ml{r,6) 



B±(B 2 - AC) 



1/2 



(17) 



where M\ denotes super- Alfvenic solution (A/ 2 > a) and 
M 2 denotes sub- Alfvenic solution (M 2 < a). The dis- 
criminant of the quadratic equation is denoted by 



B — AC 



■0 1 



pIP 2 



(GpE + GtLfie 2 



a 



(18) 



Although the locations of the Alfven point are given 
by the flow parameters L/E and Qf for each magnetic 
field line \& = constant, the distribution of the Alfven sur- 
faces in a black hole magnetosphere is obtained when the 
magnetic field distribution ^ = 9) and the boundary 
conditions at the plasma source E = E(^f), L = Lfif) 
and Qp = flp^) are specified. The plasma source (the 
plasma injection point) is located within the sub- Alfvenic 
region, where the poloidal flow velocity is less than the 
Alfven wave speed, and the neighborhood of the event 
horizon and the distant region from the plasma source 
are in the super-Alfvenic region, where the poloidal flow 
velocity is greater than the Alfven wave speed. The ac- 
cretion/wind solution is denoted by M 2 = M 2 in the sub- 
Alfvcnic region and by M 2 = M? in the super-Alfvenic 
region. Both branches of the solutions always connect 
smoothly at the Alfven point, where B 2 — AC — 0; that 
is, (M£)a = (M^a = a a- Furthermore, at the light sur- 



four-velocity of a ZAMO seen by a distant observer (see faces, we see that (M_)l — 0, while (MV)l = 2e /Al 
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is finite except the case that At, = accidentally. The 
label "L" indicates quantities at the light surfaces. 

If the coefficient A becomes zero at some location for 
a given magnetic field line specified by £ 2 = £ 2 (r; , 3>), 
the Mach number Mi of the super-Alfvenic flow solu- 
tion diverges there, while M 2 has a finite value in the 
sub-Alfvenic region. This means that such a magnetic 
field line under a given flow's parameter set is unphys- 
ical in the super-Alfvenic region. In fact, from equa- 
tions ([5]) and ® both the toroidal magnetic field B$ and 
the poloidal magnetic field B p should vanish at the lo- 
cation of A = 0, where M 2 — M 2 — > oo. (Note that, 
although {B$)a=§ = {B p )a=q — 0, the ratio (3 must have 
a finite value of /3 2 = l/z at A = 0.) Thus, such a so- 
lution (where A = in the super-Alfvenic region) can 
not be accepted as a solution of trans-magnetosonic ac- 
cretion/wind. To obtain a physically acceptable MHD 
flow, the magnetic field configuration characterized by 
such a function £ = £(r, 6) should be modified to avoid 
the appearance of A = location, or another set of flow 
parameters should be selected at the plasma source. 

For the inflow (or outflow) streaming toward the black 
hole (or distant regions), if A = on the sub- Alfvenic 
flow, the solution has a finite Mach number across this 
location. However, if A = in the super-Alfvenic flow re- 
gion, the Mach number diverges there; the super- Alfvenic 
flow is only available in the A > region. Thus, A > 
in the super- Alfvenic region is an necessary condition 
for trans-magnetosonic accretion/ wind solutions. On the 
other hand, the A < region is the forbidden region 
on the super-Alfvenic magnetosonic flow solution, where 
M 2 < is obtained. 

For the cold MHD flow, we define two characteristic 
Alfven Mach numbers related to the Alfven and fast mag- 
netosonic wave speeds (see Q) 



Mi w M) 



a 
a ■ 



(19) 
(20) 



The locations of Ml 



M| w and Mj = M| M indicate 



the Alfven point and the fast magnetosonic point, respec- 
tively. Just on the event horizon, we obtain Mp M = 0. 
Figure [1] shows the ingoing and outgoing flow solutions; 
both solutions are started from the separation surface, 
which separates the gravitational force dominated re- 
gion and the centrifugal force dominated region (see [3] ) , 
with zero poloidal velocity. These solutions cross the 
M 2 = M| w and M 
(Fig. QJ,), and the \u 
and \u r 



Mp M curves in this order 



''AW 



V / E/A(i?p/47r/ m )a 



= u r FM = v /S/A(B 2 /47r/in)(a + /3 2 ) curves 
(Fig. [T|d). Because n = constant through the flow so- 
lution, the number density n diverges at the injection 
point, where u r = 0, so that our definite radial Alfven 
wave speed and the fast magnetosonic wave speed become 
zero, (uAw)inj = ( u F]vi)inj = 0, a t the injection point. 

Note that the function £ (r, 9) gives the distribution of 
the cross section of the magnetic flux tube in the poloidal 
plane. To plot the flow solution as shown in Figure [U we 
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FIG. 1: Transfast magnetosonic inflow and outflow solutions 
(thick bold curves) started from the separation surface with 
zero velocity (labeled "inj"). (a) The square of the Alfven 
Mach number M 2 vs. radius r/m and (b) the radial 4- velocity 
of the flow u r vs. radius r/m. The solution is plotted on 
the equatorial plane (8 — n/2) in Schwarzschild geometry 
(a = 0), where the magnetic field configuration is assumed 
with £(r, 9) — 1.0. The flow parameters are given as E — 1.15, 
mQ.F = 0.7fi m ax, LQ.F = 0.6, where fi max is the maximum 
value of Q.f, and is given as the double roots of a — 0. (If 
Qf > n max , no light surfaces exist along a magnetic field 
line considered.) The broken curve shows the M 2 = M^ w 
curve (left) and the \u r \ = «a W curve (right), and the dotted 
curve shows the M 2 = Mp M curve (left) and the |tt r j = iip M 
curve (right). The crossing points of these curves with the 
flow solution labeled by "A" and "F" are the Alfven and fast 
magnetosonic points, respectively. The thin curves started 
from the light surface and approached to the A — line are 
unphysical solutions of the quadratic equation. 



need to specify the stream line of the flow; in Figure [TJ 
the flow streaming along the equatorial plane is plotted, 
but the magnetic flux tube determined by the function £ 
is not a conical shape. Figure [2] shows the ingoing flow 
solutions for rapidly rotating Kerr black hole cases. Each 
inflow enters the event horizon with a finite Mach number 
(Fig. [5^) or breaks at the A = location (Fig.J^b), where 
the Mach number of the trans-fast MHD inflow diverges. 
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FIG. 2: The square of the Alfven Mach number M vs. radius 
r/m for transfast magnetosonic inflow solutions (thick bold 
curve) on the equatorial plane (9 = n/2) in Kerr geometry 
(a = 0.8m). The flow parameters are given as E — 1.0, 
mfl F = 0.5fi max , LQ, F = -0.5 (left) and m F = 0.09 (right). 
The magnetic field configuration is assumed by equation (|A2|) 
discussed in SjA] The right-hand figure shows an unphysical 
solution, where the Alfven Mach number M 2 diverges at the 
location of A — 0. The black spots (right) indicate M 2 — oo. 
The dotted curve shows the M 2 = M AW curve, and the solid 
curve shows the M 2 = M| w curve. The crossing points of 
these curves with the flow solution labeled by "A" and "F" are 
the Alfven and fast magnetosonic points, respectively. The 
thin curves started from the light surface and approached to 
the A — line are unphysical as inflow solutions. 



In Figures [T] and [21 the functions £(r, 9) are considered as 
simple models that satisfied conditions specified at some 
characteristic radii (the details of the conditions will be 
presented in Appendix [A"]l . The general properties of 
trans-magnetosonic flows discussed here do not depend 
on detail of the functional form of £ (r, 9) . 

The poloidal velocity of a cold MHD flow can be ex- 
pressed in terms of M 2 and /3 2 (or £ 2 ) with the conserved 
quantities as 



When we consider a solution for accretion onto a black 
hole, we apply M 2 = M 2 . in the sub-Alfvenic region of 
r > ta and M 2 = Ml in the super- Alfvenic region of 
r < ta (see Fig. [Tb). By using the poloidal equation ©, 
equation (|2"Tj) can be reduced to 



(22) 



a + f3 2 



which corresponds to a physical trans- Alfvenic accretion 
solution. From the requirement of u 2 > 0, we find the 

minimum energy Eo(r, 9) at each place between the inner 
and outer light surfaces. The energy E should be taken as 
E > E (r; *) along the flow (between the plasma source 
and the event horizon). Note that, at the location of 
A = in the super- Alfvenic region (if such a point ex- 
ists), the poloidal velocity has a finite value, although the 
Alfven Mach number diverges. However, the poloidal and 
toroidal components of the magnetic field vanishes there, 
while (3 has a finite non-zero value, as mentioned in the 
first half of this section. 



IV. RESTRICTIONS ON MHD FLOWS 
A. Trans- Alfvenic Flow 

Along a magnetic field line *ff(r, 9) = constant, the lo- 
cation of the Alfven point (ta, $a) is given by M 2 = a 
and 



m F = Y A 

with the definition of a function 

Y = -G^ F /G t , 



(23) 



(24) 



u 2 p = B 2 p M 4 



P 2 (Gj,E + G t L) 2 M 4 
P l(M 2 -a) 2 



(21) 



where L = L/E. Although Lilp is a function of here- 
after, we can regard L£lp as a function of ta along the 
magnetic field line considered; note that 9a — 6*a(^a; 
In the black hole magnetosphere, two surfaces of the 
Alfven points (i.e., the Alfven surfaces) for inflow and 
outflow are distributed between the inner and outer light 
surfaces given by the relation a — 0. For an ingoing 
MHD flow, the region between the Alfven surface and 
the event horizon is the super- Alfvenic region, while the 
region between the plasma source and the Alfven surface 
is the sub-Alfvenic region. 

Now, we consider magnetic field line connected to a 
black hole with a certain value of Of- Figures^ and [5)3 
show the locations of the Alfven radii for a given LSlp 
value. When the magnetic field line rotates faster than 
the black hole (i.e., a slowly rotating black hole case), < 
^>h < S1f, whose state is named "type I" in Takahashi 
et al.Q, the condition for LQp for the existence of the 
Alfven point in the flow solution is (Lf^min < Lflp < 1 
(see Fig[3Ji), where (LQp) m i n is the minimum value of 
L£! p for the Alfven point to appear on the flow, and it is 
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given by dY A /dr A = 0. The MHD flow with E > and 
L > only is obtained. 

For the type I case, it is possible to select the inner 
or outer Alfven point in an accretion solution, and then 
trans- Alfvenic accretion solutions with two types are ac- 
ceptable (see [ill]); that is, "magneto-like" and "hydro- 
like" accretion solutions. The magneto- like accretion so- 
lution passes through the inner Alfven point and the in- 
ner fast magnetosonic point, and the hydro- like accretion 
solution passes through the outer Alfven point and the 
middle fast magnetosonic point. Note that, in the latter 
type solution case, the acceptable range of Lflp is mod- 
ified to (LfifOmin < Lflp < (Lto F ) max , where (Lfl F ) max 
is the maximum value of Lfl p to avoid the forbidden re- 
gion discussed in 0] between the outer Alfven point and 
the event horizon, and is given by dk^/drx — with 
&a = 0. For a counterrotating case afl F < (named 
"type III"), we see the similar restrictions on Lflp as 
type I, although we obtain L < and E > flows. 

Figure [3Jd shows the rapidly rotating black hole case 
of < fl F < loh (named "type II"). For a certain Lflp 
value, one Alfven point is obtained on the flow solution. 
In this case, the negative energy accretion with E < 
and L < is possible when Lflp > 1. The flow with 



E > and L < is realized when Lflp < 0. For < 
Lflp < (Lfl F ) meix , the inflow with E > and L > is 
realized. More detail discussions about the restriction by 
the Alfven points are presented by Takahashi et al.Q- 

B. Super- Alfvenic Flow 

After passing through the Alfven point, the ingo- 
ing/outgoing flow solution passes through the fast mag- 
netosonic point automatically. However, if the location 
of A(r; VP) = appears on the flow solution consid- 
ered in the super- Alfvenic region, the Mach number of 
this super-fast magnetosonic flow solution diverges there; 
that is, no physical MHD accretion/ wind solution exists. 
To obtain a physical MHD accretion/ wind solution, no 
A = surfaces in the super- Alfvenic region must be re- 
quired. Although the value of A depends on the flow 
parameters flp, E and L under a given function £(r, 9) 
and the hole's spin a, we will find restrictions on these 
combinations of parameters to avoid the appearance of 
the A = location. 

The condition A > for a super- Alfvenic flow can be 
reduced to 



A(r;#) 



-G t (l + Y + X)E 2 
(G^flp) 2 



m F - (m F ) + m F - {Lfi F ) 



> 



(25) 



with 



(Lflp^ = 



Y 



1 + Y + X 



(1+X)±e 



1 + {1-Y)X-(1 + Y + X)^- 

E 2 . 



1/2- 



(26) 



where X = g^G t (l - £, 2 )/Pw and e = ^f(^_f - 
u))/\n F (Sl F -lu)\. When the function £ 2 = £ 2 (r;1<) is 
specified as a given magnetic field structure with the field 
aligned flow parameters Q F and E, the value of (iSl f) ± 
is determined along the stream line. In Figures [3J: and 
|3Ji we show the Lflp = (if2j?) curves for the slowly 
and rapidly rotating black hole cases, respectively. The 
models of £ 2 (r; \&) are discussed in Appendix |A1 

To obtain a trans-fast magnetosonic accreting flow 
onto a black hole, the LQ F — constant line must not cross 
the LQ F = (L£lp) ± curves between the Alfven point and 
the event horizon on the Lflp - r plane. The condition 
A > for the super-Alfvenic flow solution is satisfied 
when 

(i) (Lflp)- < Lflp < {Ln F )+ 

in the G t (l + Y + X) > region, and 

(ii) Lflp < {Lflp) + or Lfl F > (Lflp)- 



in the G t (l + Y + X) < region. The value of (LSl F )- 
diverges at the location of A = — (1 + Y), while (Ln F ) + 
has a non-zero finite value there. Thus, the value of Lflp 
(or the location of the Alfven point) is restricted. Ex- 
amples of this restriction on Lflp for the A > regions 
are also shown in Figures [3J: and [3J1 for Type I and II 
cases, respectively. For accretion onto a slowly rotating 
black hole (type I; see Fig. [3);), we see that the case (i) 
is applied everywhere. For type I and III, the location 
of A = — (1 + Y) may exist, although it depends on 
a model of £ 2 ; that is, the case (ii) may arise. In this 
case, however, there is no additional restriction in the 
(LQ F ) m i n < Lflp < 1 range. On the other hand, the 
case (ii) arises near the event horizon for a rapidly rotat- 
ing black hole case (type II; see Fig. [3ji ). Note that for 
type II the case (i) is also possible just outside the event 
horizon, although it depends on £ (to be discussed in 
§ IV A[) . In this case, the negative energy inflow solutions 
only are obtained; the positive energy inflow solutions 
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FIG. 3: (TOP) The possible range of LQf for (a) type I and (b) type II cases. The vertical axis is LQf, while the horizontal axis 
is the Alfven radius. Curves LQf = Y(ta) (marked by [1a]) and k(rA) = (marked by [&a]) are plotted. Two crossing points 
between the LQf = Y(rA) curve and LQf = constant line in the type I case (left) show the inner and outer Alfven radii. On 
the other hand, one Alfven radius exists in the type II case (right). The hatched region bounded by &a = shows the forbidden 
region for the parameter. For type I case the range of (Z/f2i?)min < LQf < 1 is acceptable as a trans-Alfvenic solution, while 
for type II case the ranges of LQf < (LQp)ma,x and 1 < LQf are acceptable. (BOTTOM) Curves LQf = (LQf)^ (marked by 
[+] and [— ]) and LQf = Y (marked by [Y]) for (c) type I and (d) type II cases. Some accretion solutions with certain LQf 
values are shown by horizontal arrows, where the crossing point with the LQf = Y curve (i.e., Y — Ya) indicates the Alfven 
point (marked by the fill square). When the solution crosses the LQf ~ (LQf^ curve, where A = 0, the Mach number of the 
solution diverges there; that is, no physical accretion solution exits (the horizontal arrow marked by the cross). The shaded 
regions show the A < regions. Note that the k > region is enclosed in the A < region. In these schematic diagrams, we 
set up the field configuration £ 2 by equation (|A1|I for type I and equation (|A2|) for type II near the equator as typical models. 



are forbidden because the A < region appears on the 
way to the horizon. 

In the last of this section, we should mention that there 
is the innermost limit of the inner Alfven point (labeled 
by "A*" in Figs. [3]: and[3ji) under the given parameters 
f2i? and E. This limit gives the restriction on the Ltlp 
value, which is the maximum value for type I case and the 
minimum value for type II case. The details are discussed 
in Appendix [Bl 



V. BOUNDARY CONDITIONS AT THE 
EVENT HORIZON 

A. Inclination of Magnetic Field Lines 

Without the trans-field equation and the informations 
for the plasma sources (as the initial conditions for MHD 
flows), we can not obtain the concrete function of £(r, 9) 
in the black hole magnetosphere considered. Neverthe- 
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less, just on the event horizon, we have already known 
that = 1 (or f3 2 H = —an), which is the boundary con- 
dition at the horizon. However, just outside the event 
horizon, the magnetic field configuration depends on the 
plasma inertia effects. Here, we discuss the function 
£ 2 (r, 8) near the event horizon. One may expect that the 
signature of the function Ah{9) = A(rjj, 0) is at least de- 
termined for any field aligned flow parameter sets. How- 
ever, at the event horizon, the value of the function X 
in equation (f2"5| can not be specified by only the con- 
dition £ H = 1. Now, we will expand the function £ 2 as 
£ 2 = l+x( r / r H — 1), where x = r H(dr&H represents the 
magnetic field configuration near the horizon. To obtain 
physical MHD accretion solutions satisfying the condi- 
tion Ah > 0, we should discuss the restrictions on the 
value of x and on the allowable ranges of the field aligned 
flow parameters (e.g., E, L and Hp), which should be 
consistent with the boundary condition at the plasma 
sources. Note that the function \ can be expressed as 
the differentials of B 2 and B 2 ^ at the event horizon. 

The restrictions on the function specify how the 
inclination angles change (having poloidal and toroidal 
components) going away from the horizon. In other 
words, the x value represents the deviation of the inclina- 
tion (or pitch) angle at the horizon from the £ 2 (r, 9) = 1 



model. To consider this situation, we also expand the 
magnetic field components as B p = B P H+SB P H{r/rH — l) 
and Bcj, = B^h + SB^h^/^h — 1)- Then, we obtain 

(£V = — \ 1 + 2 ( S *Z- S **\ (il-iYI 

\Bj (-a H ) [ V Bp B<f, ) H \r H )_ ' 

(27) 

When we modify the magnetic field configuration 
(SB p /B p )h > and/or (SB^/B^h < 0, we see that the 
value of x increases. Although the term B p means the di- 
vergence of the magnetic flux tube to the (+r)-direction, 
when the ratio of the divergence (SB p /B p )h for a certain 
model is larger than that of the £ 2 = 1 (x = 0) model, 
we see that (SB p /B p )h > 0. Furthermore, if the bending 
angle to the toroidal direction is smaller than that in the 
force-free case (SB^—O along a magnetic field line), we 
see that (SB^/B^h < 0. 

If the location of A = appears between the Alfven 
point and the event horizon (i.e., the super- Alfvenic re- 
gion of the flow) , such a solution is unphysical as an ac- 
cretion solution as mentioned in § IIVB1 At the event 
horizon, to accrete onto the black hole, the condition 
Ah{x) > must be required. The boundaries of the 
Ah (x) > region on the x~LClp plane are obtained by 
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ftp + lo h + luhX h 



(1 + X h )± e H 



0| 



9u 
E 2 



Ll>h 



i - 4 1 v« 

E 2 



1/2' 



(28) 



where en = 0,p(Clp — ujh)/\91f(91f — ^h)\, Xh = 
X(r H ) = Ux and H = 2m 2 r H Gl/[{r H - ro)E*]. Fig- 
ure[Hshows the relation between (LQ,f)h an d X f° r type I 
and II cases. There is a maximum x value (= Xmax) for 
existing A > region in both type I and II cases. 

The value of {LVLf) h becomes zero at x = XOj an d 
diverges at x = Xoc- Although these characteristic x's 
depend on the flow parameters, the details are discussed 
in Appendix [C] Furthermore, the relations between the 
X value and the acceptable ranges of LQ p are also sum- 
marized there. 



B. Black Hole Spin-Up/Down via MHD Accretion 

In § IIVBI and § IV Al (see also Appendix [C]), we have 
discussed the restriction on L/E for the MHD accre- 
tion onto a black hole to avoid the A — surface on 
the flow solution, by considering the conditions at both 
the event horizon and the Alfven point. Here, we will 
discuss the increase or decrease of the angular velocity 
(spin- up/down) of a rotating black hole by applying this 
restriction on L/E to equation Hereafter, we treat 



the case of a > and Hp > mainly (i.e., type I and II 
cases). 

Although the accreting gas carries the mass and an- 
gular momentum into the black hole, both the cases 
Sloh/ojh > (spin-up) and Sujh/^h < (spin-down) 
are possible for the positive energy (E > 0) MHD in- 
flows. That is, the range of L9lp > (L£lp) ul gives the 
hole's spin-up and the range of Lflp < (L£If) u gives the 
spin-down, where 

(iO F ) w = (l + r H /m)an F ; (29) 

in Figure [4] we show the Lflp = (LD,p) u line (where 
8u>h = 0). [Similarly, with respect to the spin pa- 
rameter a, we see that 5(a/m) > (or S(a/m) < 0) 
for positive energy MHD inflows with Lilp > 2a9lp 
(or Ltlp < 2ailp).] Furthermore, by considering the 
restrictions from the Alfven point and the event hori- 
zon, we obtain that increase of the angular velocity of 
the black hole 5loh > by MHD inflows is realized 
when (LQp)^ < LQp < (ZOj?)a* for type I case, and 
(LQf) u < LQp < (LQp) H for type II case. Note that, 
even if positive angular momentum L > of MHD flow 
falls into the black hole, the angular velocity of the black 
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FIG. 4: The horizon's boundary condition on LQf for (a) type I of a — 0.3m and (b) type II of a — 0.9m. The bold thick 
curves denote the LQf = (LQf) h {x) curves (marked [+] and [— ] respectively) with E 2 = 1.0, Qf = 0.5f2 m ax = 0.09623/m 
and 9 = n/2. The sheaded region corresponds to the Ah < region, which is forbidden as a black hole accretion solution. The 
hatched regions in (a) show the LQf > 1 and LQf < regions (unphysical). The LQf = Yh line is denoted by the dotted 
line and labeled by [SSh]- The horizontal broken line marked by [Suih] shows the LQf = (LQf)u, line. Although, in (b), it 
seems (LQf)^ < (LQf)u < (LQf)h in this plot, we can find the case (LQf)lu < (LQf) h < Yh for larger (— x) values. For 
the range (LQf ) w < LQf < (LQf)~^, the spin-up is caused for type II. 



hole can be decreased; that is, 6u>h < is obtained for 
< LQ F < (LQ F ) U . 

For a slowly rotating black hole of type I case, equato- 
rial (low-latitude) inflows onto the black hole contribute 
to the spin-up, but spin-down effects can be also realized 
by MHD plasma accreting from the polar (higher lati- 
tude) region of the event horizon. The value of (LQp)min 
can be specified by the latitude of the Alfven point, where 
ta(0a) = = r A ut , so that we define the critical angle 
of 9a [= (#a)ci] by (LQ F ) min = (LQf)u>, except for a 
Schwarzschild black hole, where (LQp)^ becomes zero. 
When 9a > (0A.)cr along a flow considered, the rela- 
tion of LQp = (LfifOmin > (LQf)u is satisfied and the 
MHD inflows only contribute to the spin-up of the black 
hole. On the other hand, when 9a < (#A.)cr, the situa- 
tion (ZfifOmin < LQp < (L£If)u (the black hole's spin- 
down) is also possible; note that, even if 9a < (#A)cr, the 
MHD inflows in the range of (LQ F )^ < LQp < (LQp) a* 
can lead the black hole into the spin-up. For larger 
values of a or for smaller values of Op-, such a spin- 
down range specified by the critical value angle (6>A) C r 
expands toward lower latitude regions, while the value 
of (LVLf)ui increases; in the a — > M limit, we see that 
(LQf)u — > Yh = VIf/<jJh (> !)• Thus, the spin-down ef- 
fect dominates for a rapidly rotating black hole of type I; 
that is, even if Qf > &H> the MHD inflows spin down 
the hole's rotation. 

For a rapidly rotating black hole of type II case, both 
L > and L < inflows are possible. The inflow 
with L > decreases the spin of the black hole when 
< LQp < {LQf)u- On the other hand, the L > 



inflows increase the hole's spin when (LQf)ui < LQf < 
(LQp) H - However, after the black hole spin-up, the value 
of (LQf)ui increases, so that the parameter range of LQp 
for the spin- up decreases. Note that, in the a — * m limit, 
any ingoing flows with L > cannot increase the an- 
gular velocity and the spin parameter of the black hole 
no more, because both (LQ F ) U and 2clQf become Yh, 
where (LQ F ) H < Y H < (LQ F ) H ; that is the A H > 
region. On the other hand, in type II case, the negative 
angular momentum MHD inflows are also possible. In 
such a case, we always obtain 5u>h < and 5(a/m) < 
as expected. Furthermore, if we consider the negative 
energy (E < 0) MHD inflows [7[, which always carry the 
negative angular momentum (L < 0) , we always see that 
Sloh/o^h < and 6(a/m) < 0. 

Figure shows the summary of the hole's spin- 
up/down by MHD inflows, which also includes the con- 
tribution by the magnetic field parameter \- First, 
in the case of type I (see Fig. [S^,), the spin- up of 
8loh/uh > is obtained when the value of LQp is in the 
range of MAX[(LQ F ) min , (LQ F )u,, (LQf) h ] < LQ F < 
(LQp) a* < 1- The spin-down of 5ojh/^h < is obtained 
when MAX[(LQ F ) min: (LQp ) H ] < LQp < {LQ F ) U , 
which is only realized for inflows with 9a < (^a)ci- 
Next, for type II case (see Fig. [5}d), we obtain the spin- 
up of Suih/^h > when (LQp)^ < LQp < (LQp) H 
for x < Xwi where x^ i s defined by the condition of 
(LQp) H = (LQp)u- The spin-down of Sujh/^h < 
is obtained when (1) LQp < MIN[(LQ F ) U , (LQf ) h } or 
(LQ f )a* < LQp for x < Xoo and (2) (LQ f )a* < LQp < 
(LQf)h f° r Xoo < X < Xmax- When we consider a coun- 
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FIG. 5: Summary of the LQf ranges for the spin-up/down of a black hole in (a) type I case and (b) type II case. The value 
of LQp for acceptable MHD accretion is restricted by at least the boundary condition at the event horizon (Ah > 0) and the 
condition of the Alfven point for existing in the A > region. The shaded forbidden region is the same as that of Fig. H The 
labels [A*], [Su>h] and [SSh] show the LQf = (LQf)a* line, the LQf = (LQf) u line and the LQf = Yh line, respectively. The 
hatched region is forbidden by the condition at the Alfven point. For type I case, the range of (LQf)hAvl < LQf < (LQf)a* is 
acceptable, and the MHD inflows with 9a > (0\)cr always cause the hole's spin-up. 



terrotating magnetosphere of Sip < (type III), where 
a > is considered, we have the relation LSlp > > 
2aSlp > (LSlp)^. In this type III case, both the angu- 
lar momentum and the spin parameter of the black hole 
decrease; that is, Sloh/^h < and 5(a/m) < 0. 

From the definition of the surface area of the black 
hole Sh = 47r(r| f + a 2 ), we obtain the relation SSh oc 
E — ujhL. Here, we should remember that, in the case of 
type II, the acceptable range of LSlp for positive energy 
(E > 0) MHD inflows is given as LSlp < (LSlp) H < Y H , 
while LSI p < 1 < Yh in the case of type I and LSI p > > 
Yjj in the case of type III. So, for E > flows (in all cases 
of types I, II and III), we can show SSh > 0, which means 
the increasing surface area of a black hole. On the other 
hand, for E < flows (that can be possible in the case 
of type II), the relation of Yh < LSI is always satisfied, 
and then the relation of SSh > is also satisfied. Thus, 
the area low SSh > of the black hole is confirmed for 
stationary and axisymmetric ideal MHD accretion flows. 



C. Secular Evolution of the Black Hole Spin 

In § IV B[ we mentioned that the MHD inflow prevents 
the evolution to the maximally rotating black hole. This 
is because Sujh < and S(a/m) < for the extreme 
rotating black hole limit. On the other hand, we know 
that Sloh > and S(a/m) > for the non-rotating black 
hole case. Now, we will discuss a secular evolution of 
the black hole spin by MHD accretion. Here, we assume 
that the values of L(^f), E^) and SIf^) keep constants 
for a magnetic flux tube considered (a <iI/=constant line) 
during the evolution. Furthermore, we do not discuss the 
plausible configuration of the magnetic field as a global 



solution of the black hole magnetosphere. 

First, we consider the case of type I (0 < ujh < Sip). 
When the magnetized plasma of (LSlp)^ < LSlp < 
(LSIp)a* falls into the black hole, the spin of the black 
hole will increase in secular time-scale. Then, the value of 
[LSlp)^ also increases, and the state of LSlp = (LSlp) u , 
where Sloh = 0, will be realized sometime. In this state, 
the evolution of the black hole's spin is terminated. On 
the other hand, when (LSlp) min < LSlp < (LSlp)^, the 
black hole's spin decreases, and the value of (LSlp)^ 
is also decreases. The final state will be settled to 
LSlp = (LSlp)^. Although the spin value a terminated 
by accreting MHD plasma is given by (LSlp)^ = LSlp, 
the values of L, E and Sip should be specified by a magne- 
tized accretion disk model. Note that, during the spin-up 
stage, the situation of luh —> Sip may be realized before 
the LSlp = (LSlp)^ state. In this case, the spin value a 
terminated is specified by Sip. 

Next, we consider the type II case (0 < Sip < u>h)- 
When magnetized plasma with (LSlp)^ < LSlp < 
(LSlp) H , the black hole will spin up in secular time scale. 
In this case, the value of (LSlp)^ also increases and the 
situation of LSlp = (LSlp)^ would be realized. When 
< LSlp < (LSlp)^, the black hole spins down in spite 
of the positive angular momentum inflows. With decreas- 
ing the hole's spin, the value of (LSlp)^ also decreases, 
and the state of LSlp = (LSlp)^ would be realized. One 
may expect that, by decreasing the angular velocity of 
the black hole, the state of luh — > Sip is achieved before 
reaching the state of LSlp = (LSlp)u- However, the state 
of u>h —> Sip is only possible for the L < inflows (see 
Appendix [Cjl . So that, the inflows with L > are termi- 
nated to the state of LSlp = (LSlp) u . In type II case, the 
negative angular momentum inflows are possible when 
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LVlp < or LVLp > {LVLf)a*- The rotating black hole 
spins down, and reaches the state loh — ► Sip. For the 
magnetically dominated accretion flows, where E < 
and L < 0, such a state would be only realized after a 
long time spin evolution. Although the state u>h ~~ * ^f 
can be achieved in both type I and II, the transition 
from type I to type II or the reverse is impossible. This 
is because the asymptotic evolution ujh — > &f is possi- 
ble for L > inflows of type I or for L < inflows of 
type II. The transition between type I and II contradicts 
the model's assumption that the field-aligned parameters 
do not change during the hole's evolution. 

In this section, we have discussed about the spin- 
up/down effects along one magnetic flux tube. However, 
in actual black hole magnetosphere models, the above ar- 
guments could be considered for each magnetic flux, and 
the effects on the black hole spin should be integrated 
over the whole magnetic fluxes across the event hori- 
zon. In type I case, although the equatorial inflow would 
dominate in the black hole accretion and such an inflow 
contributes to the spin-up, but the spin-down by the in- 
flow in the polar region may not be negligible in a black 
hole magnetosphere. We must consider the spin-down 
effect carefully in such a case. In the polar region, the 
state (Lf^min < {LQp) u is easily obtained except for 
a non-rotating black hole case. Then, the ingoing flows 
can carry the energy into the black hole with less angu- 
lar momentum; this means the spin-down of the black 
hole. When disk's gases fall into the black hole along the 
disk-black hole connected dipole-like magnetic field lines 
[l2l [III], such a situation would be possible. Thus, the 
configuration of the magnetic field would be very impor- 
tant in astrophysical situations. To estimate the actual 
spin-up/down of the black hole, we need to integrate the 
spin-up/down effects on each magnetic flux tube from the 
pole to the equator. To do this, it is necessary to con- 
struct a reasonable model of magnetized accretion disk 
for distributions E = E(0 H ), L = L(0 H ), M F = ^f(Oh) 
and rj = i~i{6h). Furthermore, the function \ = x{@h) 
over the event horizon should be obtained from the stud- 
ies of the GS equation around the event horizon. Then, 
we can estimate the evolution of the rotating black hole. 
However, this problem remains as our future task. 



VI. CONCLUDING REMARKS 

We have discussed about stationary and axisymmetric 
ideal MHD inflows onto a black hole. We have shown 
that a non-singular distribution of /3(r, #) gives ingoing 
and/or outgoing trans-magnetosonic flow solutions with- 
out critical behavior. Such solutions automatically sat- 
isfy the regularity conditions at the magnetosonic point, 
and passes through the magnetosonic point. In a cold 
flow case, we only solve the quadratic equation for M 2 
to obtain trans-fast magnetosonic flow solutions. The in- 
going flow passes through the inner or middle fast mag- 
netosonic point. When we consider a hot MHD flow, we 



can also treat the flow solution as a polynomial of high 
degree in M 2 without the regularity condition at the fast 
and slow magnetosonic points. 

In this paper, to discuss cold MHD inflows onto a black 
hole, we have not specified the function /3 (or £) without 
the GS equation. However, we have found the restric- 
tions on the field-aligned flow parameters under a given 
field geometry. The ranges of possible LVlp values for 
black hole accretion are restricted by the condition at the 
Alfvcn point as discussed in § IIV Al and Appendix [B] and 
by the condition A(r; ^) > along the magnetic stream 
function that is discussed in § IIVB1 where the function A 
depends on the MHD flow energy E, the field geometry 
£ (or /3), and the angular momentum related parame- 
ter Lfip. Although the value of Lflp is related to the 
energy and angular momentum of the ingoing flow, the 
field geometry around the event horizon restricts these 
values. Furthermore, we have discussed the inclination 
X of magnetic field lines at the event horizon. When 
we try to solve the magnetic field distribution from the 
plasma source to the event horizon in a black hole mag- 
netosphere, we should be aware of this boundary condi- 
tion. Note that the restrictions and conditions discussed 
in this paper are model independent as we are focusing on 
stationary and axisymmetric black hole ideal MHD ac- 
cretion. In future studies of black hole magnetospheres 
(including numerical studies), our approach will be help- 
ful to check the magnetic field structure near the event 
horizon. 

With the restriction on L/E for MHD accretion, we 
have discussed the secular evolution of black hole's rota- 
tion for certain magnetic flux-tubes. Although we have 
applied the boundary condition at the event horizon, our 
approach by the £-model guarantees that the ingoing 
MHD flow onto a black hole is trans-fast magnetosonic. 
Then, we find that there are two asymptotic states for 
the spin evolution. One is the state that (i) the angular 
velocity of the black hole approaches to that of the mag- 
netic field line, loh — ► fij?, and the other is the state that 
(ii) the spin-up due to angular momentum influx and the 
spin-down due to ingoing mass influx of MHD flows are 
canceled, where LQp = (LQp)^ is achieved. For type I 
case the asymptotic state of Sujh — is achieved when 
Lflp = (LQp)^ (see Fig. [5k), although the state (i) of 
luh — ► Of may be realized before the state (ii). In the 
case of type II (see Fig. [SJd), when L > the asymp- 
totic state is given by Lflp = (LQp)^, while the spin 
terminates as u>h — > ^f when L < 0. 

When we discuss the actual asymptotic state in as- 
trophysical situations, we should integrate the magnetic 
flux-tubes over the event horizon. Although we can spec- 
ulate on the final stage that the angular velocity of the 
black hole will be settled to a typical value of that of the 
whole magnetic field lines connected to the event horizon, 
we need to construct a reasonable MHD accretion model 
under a certain magnetic field configuration to specify the 
functions L(^>) and E(^). The basic idea of the black 
hole-disk connection of magnetic field lines has been sug- 
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gested by some authors (e.g., 0, E3, El, El HJ EH). 
The realistic features of magnetic field lines have been 
presented by [H, EE EI 0, El HI ■ In an astrophysical 
point of view, the magnetic connection between the black 
hole and the disk would play a very important role. The 
power input /output from the black hole depends on the 
ingoing MHD flow properties and the shape of the mag- 
netic field lines. The released energy of the plasma in the 
deep gravitational potential well can be carried directly 
from the distant disk surface to the horizon by the black 
hole-disk connecting magnetic field lines. Furthermore, 
when the black hole is rapidly rotating, the rotational en- 
ergy of the black hole can be carried from the black hole 
to the disk through magnetic interactions. Such energy 
and angular momentum transport from the disk to the 
black hole determines the fate of hole's rotation. 

Although in this paper a stationary magnetosphere is 



discussed, their dynamical phase is also an important 
problem. The structure of a magnetosphere around the 
black hole is also discussed in [2y, [2?], El, Eil by general 
relativistic MHD numerical simulations. Our current in- 
vestigations performed by analytical methods can help 
gain deeper insight from the results obtained by time- 
dependent general relativistic MHD simulations. 
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APPENDIX A: BENDING ANGLE OF 
MAGNETIC FIELD LINES 



We will discuss the conditions on the electric-to- 
magnetic field amplitude function £ or the bending angle 
(3 along a magnetic field line connected between the disk 
surface and the black hole. For a given function £ 2 = 
^ 2 {r; \&) or (3 2 = f3 2 (r;^i) we can easily obtain a trans- 
fast magnetosonic accretion solution M 2 = M 2 (r;^) 
by the quadratic equation as discussed in § IIHI To 
obtain a physically acceptable accretion solution, how- 
ever, we should discuss some requirements on the func- 
tion £ at some characteristic radii, where some condi- 
tions may restrict the functional form of £ 2 . To discuss 
these requirements on £ 2 , we will consider a disk - black 
hole connecting magnetic field line, where the plasma in- 
jected from the disk surface streams toward the black 
hole along a magnetic field line. First, the magnetic field 
line would almost corotate with the "footpoint" on the 
plasma source, where the magnetic field line is anchored. 
So we may expect B$ ~ ((3 2 3> 1) there when the 
toroidal surface current dominates the poloidal one, or 
f3 2 ~ 0(1) when the toroidal surface current is almost 
the same as the poloidal one. 

In the cases of type I and III, the ingoing plasma in- 
jected from the disk surface may bend the magnetic field 
toward the counterrotating direction (the trailed-shape of 
the magnetic field line: < 0) due to the plasma iner- 
tia effect. This situation is analogous to pulsar magneto- 
sphere, although the magnetic field lines extend inwardly 
to the black hole. Conversely, the plasma flow ejected to- 
ward the rotation direction (the leading-shaped: B^ > 0) 
may be also possible. In this case, however, the magnetic 
field line flips over at the "anchor point" (see also [30l|). 
which corresponds to the Alfven radius of M 2 ^ a (where 
B0 = 0). At the anchor point, we see that £ 2 = oo and 
(3 2 = 0. At the Alfven point that is the Alfven radius 
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TABLE I: Restrictions on the functions £ 2 and f3 2 at various 
characteristic radii are summarized. The radius appeared in 
accretion solution is marked by "o" for types I /III and II, 
while the mark "x" shows the absence of that radius in the 
solution. 



Characteristic Radii 


e 


(3 2 


I/III 


II 


Event Horizon (H) 


TH 


i 


— OIH 


o 


o 


Corotation point (C) 


r u 





finite 


X 


o 


Alfven Point (A) 


TA 


finite 


finite 


o 


o 


Anchor point (AN) 


TAN 


CO 





o 


X 


Plasma source (I) 


ri 


finite 


finite 


o 


o 



with the condition M 2 = a, we see B^ ^ and both £ 2 
and P 2 have finite values. 

Next, in type II case, the rotating black hole must 
bend the magnetic field toward the rotating direction 
(Bfj, > 0). In this type, no anchor point exists, but the 
"corotation point" where Qf = w exists between the in- 
ner and outer light surface, while in type I and III cases 
no corotation point exists on the way to the event hori- 
zon. The value of £ at the corotation point, where Y = 0, 
must become zero (i.e., 1 + X = 0) to avoid the A < re- 
gion near the corotation point, while f3 2 has a finite value 
there. [If not so, we obtain (LfL F ) + = (LClp)~ = at the 
corotation point, and we see the A < region between 
the plasma source and the event horizon; that is, no phys- 
ical ingoing solution exists.] Finally, at the event hori- 
zon, the boundary condition there requires that £|f = 1 
or prj = —an- These requirements on £ at characteristic 
radii are summarized in Table |TJ 

For example, for type I and III cases, we can simply 
set a model with £ 2 (r, ff) = £|j = 1 throughout the flow 
region considered for ingoing flows without the anchor 
point. (Here, we consider a situation that the plasma 
source is placed inside the anchor point.) As another ex- 
ample, we may add a deviation term from the £ 2 (r, 9) = 1 
model, which must become zero on the event horizon, as 
follows; 



e = i + c-i 



(Al) 



where C\ is a constant. 

To see the behaviors of accretion onto the black hole in 
type II case, however, we should consider the function £ 
to satisfy both the boundary condition at the event hori- 
zon and the requirement on the corotation point. Then, 
for example, we will consider the following function for 



e = 



( LU — flp 
\UJ H - ftp 



(A2) 



where C\\ is a constant. [The functions (LQp)^ seen in 
Figure [5Ji could be plotted under this distribution of £.] 
Note that the £ 2 (r, 9) = 1 inflow model in the type II case 
requires (3 = at the corotation point. This requirement 



means = that corresponds to the anchor point, but 
there is no anchor point in type II flow solutions. Thus, 
the simple £ 2 (r, 6) = 1 model can not be applied as an 
inflow solution for a rapidly rotating black hole (type II) 
case. 

For outflows (r 3> to), for example, we can make a 
model 



(A3) 



where Co is the parameter for the magnetic field geometry. 
For Co < 0, the outgoing flow reaches the distant region 
with a finite Mach number, while for Co > the flow is 
confined within a certain radius R c where the Mach num- 
ber becomes to diverge (see 13], for detail discussions). 

In a realistic situation for the black hole magneto- 
sphere, a global structure of magnetic field lines would 
not be able to be expressed by simple forms of func- 
tion £. To understand the basic properties of MHD in- 
flows/outflows, it will be helpful to use function (|A1|) or 
(|A2|) for ingoing winds and function (|A3|) for outgoing 
winds. In § IIII1 we show Figures [T] and [5] that are ac- 
cretion solutions along a given stream line with a given 
function £(r;W). However, we should note that, as a 
practical matter, we can obtain the velocity distribution 
u r (r, 9) in the poloidal plane when the function £(r, 9) is 
specified as a model. 

Now, we consider the restriction on the function £(r, 9) 
for the A > region, which is the necessary condition for 
MHD ingoing/outgoing winds, in the black hole magne- 
tosphere. If the discriminant in equation (f2l)]) becomes 
zero at some location (r, 9) in the super- Alfvenic region, 
equation pi?)) becomes 



(m F ) + = (Ln F y = {m l 



(A4) 



where 



(Ln F ) d = l - |i , (A5) 

and the corresponding £ 2 value, is expressed as 



g<p<p(QF ~ u) 2 (g tt - E 2 ) 



aE 2 



(A6) 



Under the condition of £ 2 (r; "P) < £^(r; ^) along a mag- 
netic field line \& = ^(r,9), the requirment A > on 
the super- Alfvenic solution is possible for a suitable LVi F 
value at most in the magnetosphere. If the situation 
with £ 2 (r;\I>) = £ 2 is achieved on the magnetic field 
line considered, the A = surface must appear on the 
way to the £ 2 = C 2 . radius (i.e., between the plasma 
source and the location with £ 2 = £ 2 ). Such a situa- 
tion gives an unphysical flow solution. Note that, for 
the inflow, if the £ 2 = C 2 , surface encloses the event hori- 
zon, no black hole accretion solution exists. Although 
the value of £^ diverges (£ 2 = ±oo) at surfaces with 
a/G 2 = l/E 2 [ or Y = (LQ F ) d ], the inner and outer 
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surfaces of £, 2 {r,6) = oo exist in the black hole mag- 
netosphere. Between these surfaces, the value of ^ is 
always negative and the discriminant in equation (|26l) is 
positive; that is, the condition A > is achieved every- 
where. This £ 2 < region occupies the magnetosphere at 
most between the inner and outer light surfaces. For the 
E ^> 1 flow, the £j (r, 9) — oo surfaces coincide with the 
light surfaces (given by a = Gf/E 2 ~ 0). Thus, for ac- 
cretion/wind solutions, the condition £ 2 (r; < ^(r; "I") 
must be required along a magnetic field line considered, 
except for the Q < region. 



magnetically-dominated flow limit, we see that E — > oo, 
Lf2i? — * 1 and ta* — * r™.) For the magnetically- 
dominated flow, the deviation from the force-free model 
with respect to the LVL F value is of the order of E~ 2 , 
and the location of the Alfven point also separates from 
the inner light surface in the same order [ (fA*/^i) = 
1 + \0(E~ 2 )\ ]. When we discuss some problems (e.g., the 
energy and angular momentum transport) in the MHD 
framework, the deviation from the force-free model be- 
comes important. This necessary condition for the trans- 
Alfvenic flow solutions discussed above also applies to the 
outflows. 



APPENDIX B: THE RESTRICTION ON THE 
ALFVEN POINT BY THE A > CONDITION 

Here, we consider the restriction on the Alfven point 
due to the condition A > for the MHD accretion solu- 
tion again (see § IIV Ap . Although the relation Lflp = Ya 
gives the Alfven surface, the Alfven surface must be lo- 
cated within the A > region. Then, the innermost limit 
of the inner Alfven point (marked by A*) for a given mag- 
netic flux surface is given by Ya* = (LQp)^ for type I 
and III, and Ya* = (LClp )7 !t for type II, whose relation 
can be reduced to 

(m F ) A * ee Y A * = 1 - . (Bl) 

Note that this relation is independent of the function £. 
If the Alfven point is located inside this critical loca- 
tion, the location with A = appearers in the super- 
Alfvenic region. Then, the requirement for the accept- 
able LVt p range discussed in § IIV Al should be modified 
to avoid such a situation (see, Figs [3fc and [3ji, where 
the location of r = ta* is given by the fill circle on the 
L£l p = Y curve) . That is, we should take the range with 
(LCIf)a* < Ltlp for the E < inflows in the type II case, 
where (LQ f )a* > 1- On the other hand, we should take 
the range with L£l F < (LHf)a* for the inflows (E > 0) 
in the type I and III cases, where (LSI f )a* < 1- 

Thus, with respect to Lttp for the black hole accreting 
flows, we summarize the necessary condition with 

|1-£0,|>K^1. (B2) 

If not so, A = is realized between the inner Alfven 
point and the event horizon; that is, the Mach number of 
the trans- Alfvenic solution diverges there (unphysical). 
When we consider the trans- Alfvenic ingoing flow, for 
Ta > f*A* equation (|B2[) can be rewritten as e 2 > a a, 
which gives a physical solution. However, no physical 
trans-Alfvenic flow, where u 2 (rA) < 0, is obtained for 

TA < TA*- 

When we treat a magnetically-dominated flow with 
£» 1, we should mention the deviation from the limit 
of large magnetization; i.e., the force-free case. (In the 



APPENDIX C: THE NATURE OF THE 
PARAMETER x 

In IV Al we introduced the function x that represents 
the magnetic field configuration near the event horizon. 
Although the function x(@h) should be given by solving 
the GS equation, in this paper we only treat it as a pa- 
rameter without the GS equation. Here, we discuss the 
characteristic nature of the parameter x- 

First, we consider the maximum x value (ee x ma x) for 
existing A > region in both type I and II cases, where 
we obtain 

i + mi-ui)i^M/E 2 ) 

Xmax nmF-^/^ii-gS/E 2 )' 

which gives (Lflp)^ = (Lflp)^. In the Schwarzschild 
black hole case, we have Xmax ~ (2mUp sin0) -2 , which 
is finite near the equator. When x > Xmax, no physical 
trans-fast magnetosonic accretion solution exists for any 
LQp values. On the other hand, we have that [L^tp)^ ~ 
(LQp)]j ~ Vtp/uj H = Y H when (— x) > 1. 

Second, we see that the value of (LQp)^ diverges when 

x = Xoc = -(LUH + n F )/(Hcu H ) , (C2) 

while (LQp)h is finite. For a rapidly rotating black hole 
case (type II) this situation, x = Xoo, arises in the x > 
region (see Fig.Q})). For a slowly rotating black hole case 
(type I), however, the value of Xoo is always negative, 
and then no restriction by the Ah < condition exists 
on the acceptable L£l p range (see Fig. 0J,) • Furthermore, 
we also specify the third value 

Xo = -[(r H - m)Z H /{2m 2 r H )]/E 2 , (C3) 

where the function (LS}p) H becomes zero. Although the 
function (L£lp) H has a E dependence, the dependence is 
weak for E > 1. In the limit of E = oo, we see that xo = 
0, where (Lfl F ) H = 0, and Xmax = XooOJ 2 H /{tip - 
where (LQ F ) H = 1, while Xoo is independent of E. 

In i jlVBI we introduce the restrictions on Lflp by the 
condition (|B2jl at the Alfven point. The additional re- 
striction from the boundary condition at the event hori- 
zon should be also imposed to avoid the A — surface 
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on the super- Alfvenic accretion flow. Then, the ingoing 
flow solution can be classified by the \ value. First, for 
the slowly rotating black hole case (type I) , the ranges of 
Lflp < (Lfl F )~^ and Lflp > (Lfl F )~^ are forbidden from 
the boundary condition at the event horizon. In addi- 
tion to this restriction, the range of Lfl F < (Z/£l J p) m i n 
and (Lflp) a* < Lfl F are also forbidden from the con- 
dition of the Alfven point, where (Lflp) a* < 1 for 
type I case. When x < Xmax, the acceptable range of 
Lflp is given by MAX[(Lfl F ) mm , (L^f)h\ < Lfl F < 
(Lflp) a* for the magneto-like accretion solution, while 
that is MAX[(Lfl F ) miu ,(Lfl F ) H } < Lflp < (Lfl F )+ in < 
(LD, f )a* for the hydro- like accretion solution, where 
{Lflp)min is the minimum of the function (Lfl F ) + = 
(Lfl F )+(r) (see also Fig. [3t). Note that [Lfl F ) H < 
for x < Xo, while the value of (Lfl F ) mm is always posi- 
tive. For type III case, we can see the similar behavior 
of Lflp. 

Next, for a rapidly rotating black hole case (type II), 
the value of x is related to the signature of the total 
energy E and angular momentum L of the flow. When 
X < Xoo the accretion solution in the range (Lflp)^ < 
Lflp < (Lflp)jj must be forbidden from the boundary 
condition at the event horizon, and when Xqo < X < 
Xmax the solution in the ranges Lflp < (Lflp)^ and 
Lflp > (Lflp)^ must be forbidden. Furthermore, the 
range (Lflp)~^ < Lflp < (Lflp)A*, where (L£If)a* > 1 
for accretion solutions in type II, is also forbidden by 
the existence of the A = surface between the Alfven 
point and the event horizon as mentioned in ijlVBI Then, 
we obtain the acceptable ranges as Lflp < (Lflp)~t and 
Lflp > (LQf)a* for ingoing flows with x < Xoo, while 
the acceptable range is (Lflp) a* < Lflp < (Lflp)^ for 

inflows With Xoo < X < Xmax- 

Here, we will consider three ranges of %- value for type 



II case; that is, 

(a) Xoo < X < X max • 

(b) Xo < X < Xoo, 

(c) X < Xo- 

In all cases, the negative energy inflow is acceptable. In 
the cases (b) and (c) the positive energy inflow is also 
possible, while in the case (a) energy E must be nega- 
tive. Thus, the value Xoo gives the maximum x value 
for the positive energy input to the black hole. Simi- 
larly, for the cases (a) and (b), the angular momentum 
of the inflows must be negative, while in the case (c) 
both the positive and negative inflows are available. The 
X = inflow is included in the case (b); the example 
of (Lfl F ) = (Lflp) I (r) curve shown in Figure [Sb corre- 
sponds to this case. When we consider the a — > m limit, 
where Ti. — > oo, we obtain (Lflp)f^(x) ~^ Yh — 2mfl F , 
and Xmax -> 0, Xoo — ► and xo -> 0. Then, only 
X < is acceptable, which gives the acceptable Lflp 
of Lfl F < Y H and (Lfl F )A* < Lflp. When ojh — * fl F , 
we have H = 0, Xmax -> oo, Xoo -> oo, (Lfl F ) H (x) = 1 
and (Lfl F )j { (x) — 0. The acceptable Lfl F ranges are 
Lflp < and (Lflp) a* < Lflp] only negative angular 
momentum flows are available. For the off-equatorial in- 
flow of 9 <C 1, we have H oc sin 2 9, \xoo\ oc sin -2 9 and 
Xmax oc sin~ 2 9, while the value of Xoo does not change 
drastically. Just on the pole, we find that (Lflp)^ ~ 
and that for the type II case only the inflow with L < 
is available. Thus, the ratio of the toroidal and poloidal 
magnetic fields that is related to the bending angle of the 
magnetic field line is restricted by the total energy and 
angular momentum of accretion onto a rotating black 
hole. 



